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Abstract 

In this work, we consider the canonical charmonium assignments for Y(4360) and Y(4660). 
Y(4660) is good candidate of 5^Si cc state, the possibility of Y(4360) as a cc state is studied, 
and the charmonium hybrid interpretation of Y(4360) can not be excluded completely. We evaluate 
the e^e~ leptonic widths. El transitions. Ml transitions and the open flavor strong decays of 
Y(4360) and Y(4660). Experimental tests for the charmonium assignments are suggested. 

PACS numbers: 12.39.-x,13.20.Gd, 13.25.Gv,14.40.Lb 
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I. INTRODUCTION 



In the last five years we have witnessed a revival of interest in charm spectroscopy, the B- 
factories( Babar and Belle ) and other machines have reported a large number of new states 
with hidden charm: /ic(lP)fl, ??c(2S)|2|, X(3872)3, X(3940)|4|, Y(3940)|5|, Z(3930) ^ and 
Y(4260)[3|. Some of them can be understood as cc states, while a conventional assignments 
for some are elusive( for a recent review, see, e.g. [sl ). These discoveries are enriching and 
also challenging our knowledge for the hadron spectroscopy, and the underlining theory for 
strong interactions. 

Recently the Belle collaboration has observed two charmonium-like states Y(4360) and 
Y(4660) in e+e" 7t+it-^/j{2S) via initial state radiation^. The mass of Y(4360) is 4361 ± 
9 ± 9MeV/c^ with a width of 74 ± 15 ± lOMeV/c^ and the statistical significance is of more 
than 8(T. The mass of Y(4660) is 4664±ll±5MeV/c2 with a width of 48±15±3MeV/c2 and 
statistical significance 5.8cr. Both these two structures are known to be produced in initial 
state radiation from e^e~ annihilation and hence to have J^^ = 1 . They were seen in the 
decays with products of 7r+7r~?/'(25'). It has been determined by the Belle collaboration that 

r(y(4360) ^ e+e-)i3(F(4360) 7r+ tt- if; {2 S)) = 10.4 ± 1.7 ± 1.5(11.8 ± 1.8 ± lA)eV/c'^ 
r(F(4660) ^ e+e-)i3(y(4660) ^ 7r+7r-z/;(2^)) = 3.0 ± 0.9 ± 0.3(7.6 ± 1.8 ± 0.8)eV/c2(l) 

where the numbers in the bracket are the solution II fit performed by Belle. In order to 
understand the nature of Y(4360) and Y(4660), it is worth to note that the Babar collab- 
oration has observed a broad structure Y(4325) in the process e~^e~ — > '-fisR'n'~^'n'~ip{2S) at 
4324 ± 24MeV/c2 with a width 172 ± 33MeV/c2Q. The mass of Y(4360) is close to that 
of Y(4325), the main difference between Y(4325) and Y(4360) is their widths, and it seems 
very difficult to observe these two structures simultaneously because of the large width of 
Y(4325). If both Y(4325) and Y(4360) are not experimental artifacts, they could be the 
same structure and the width difference is due to the experimental error, or they are two 
different resonances. If one assumed that they were two different structures, it would be 
very difficult to assign both of them as conventional charmoniums simultaneously. So, at 
least one of them should be exotic state. It possibly may be produced by DiD( or DsiDg 
) rescattering effect or some other mechanism, therefore it would indicate the necessity of 
refinements in the naive "quenched" qq quark models or the inclusion of additional dynam- 
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ical effects. In this paper, we assume Y(4325) and Y(4360) are exactly the same resonance 
for simphcity. Possible non-charmonium assignment of Y(4360) and the relation between 
Y(4325) and Y(4360) will be considered in detail in future work. Finally we would like to 
mention that the Belle collaboration claims the broad Y(4325) is comprised of two narrower 
peaks Y(4360) and Y(4660)|ll|. 



In order to understand the structure of Y(4360) and Y(4660), i.e. whether they are con- 
ventional charmonium states or other exotic structures, it is very necessary to first consider 
the canonical charmonium assignments and the characteristic signals. With J^^ = 1 , a 
conventional cc state is either a S-wave state or a D-wave state. There are already reasonably 



well established cc candidates for IS, 2S, ID, 3S, 2D and 4S[12l], therefore new 1 charmo- 
nium states can only belong to 3D, 4D, 5S or 6S, and a natural assignment for Y(4360) will 
be a 3 ^Di cc state, and Y(4660) as a 5 ^Si charmonium. However, this assignment has the 
problem that the mass of the Y(4360) is somewhat lower than the nonrelativistic potential 
model prediction for the 3 ^Di cc state, which will be shown in the next section. 

In this work, we study the properties of Y(4360) and Y(4660) under the hypothesis of 
Y(4360) as a 3 '^Di cc state and Y(4660) as 5 ^Si cc state. We first briefly review the nonrela- 
tivistic potential model and give its prediction for the masses of 3 "^Di, 4'^Di, 5 '^Si and 6 ^Si 
cc states. The e"'"e~ leptonic widths. El transitions. Ml transitions and open charm strong 
decays of both Y(4360) and Y(4660) are studied in section III and IV respectively. From 
these results, we suggest adequate measurements which can verify the canonical charmo- 
nium assignments and distinguish the cc structure from other non-cc possibilities. Finally 
we present our summary and some discussions. Possible charmonium hybrid assignment of 
Y(4360) and its crucial decay modes are suggested. 

II. REVIEW ON NONRELATIVISTIC POTENTIAL MODEL AND THE CHAR- 
MONIUM ASSIGNMENTS FOR Y(4360) AND Y(4660) 

The quark potential models have successfully described the charmonium spectrum, which 
generally assumes shorted-ranged color coulomb interaction and long-ranged linear scalar 
confining interaction plus spin dependent part coming from one gluon exchange and the 
confining interaction. The potental model is closely related to QCD, which can be derived 
from the QCD effective field theory Q, 3. Here we shall use the simple nonrelativistic 
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potential model proposed by T. Barnes and S.Godfrey and E.S.Swanson|15l|. the zeroth-order 
Hamiltonian is, 

Ho = — --- + br + -^5.(r)Se.S, (2) 

where Sa^{r) = {cr/y/n)^ e"^^^^ , which is a gaussian-smeared hyperfine interaction. Solu- 
tion of the Schrodinger equation with the above Hq gives our zeroth order charmonium 
wavef unctions. The splitting within the multiplets is then determined by taking the matrix 
element of the spin- dependent Hamiltonian between these zeroth-order wavef unctions. 
The spin-dependent Hamiltonian is taken from the one-gluon-exchange Breit-Fermi Hamil- 
tonian (which gives spin-orbit and tensor terms) and an inverted spin-orbit term, which 
follows from the assumption of a Lorentz scalar confining interaction. The Hsd is as follows, 

= 4 " ^' ■ 

This simple potential model consists of four parameters: the strong coupling constant 
as which is taken to be a constant for simplicity, the string tension b, the charm quark 
mass rric, and the hyperfine interaction smearing parameter a. Fitting the masses of the 
11 reasonably well established experimental charmonium states, the values of these four 
parameters are already fixed as follows: a, = 0.5461, b = 0.1425GeV^ = 1.4794GeV and 



a = 1.0946GeV 



Solving the Schrodinger equation with the zeroth-order Hamiltonian 
Hq numerically by the Mathematica program |l^ and treating the spin-dependent terms Hgd 
as mass shifts by the leading order perturbation, we obtain the masses and wavefuntions of 
the canonical cc states. The masses of 3^Di, 4'^Di, 5'^Si and 6'^Si are predicted as, 

M(3 ^Di) = 4455MeV, M(4 ^Di) = 4740MeV, 

M(5 ^Si) = 4704MeV, M(6 ^Si) = 4977MeV (4) 

Comparing with the masses of Y(4360) and Y(4660), it is natural to assign Y(4360) 
as a 3^Di and Y(4660) as a 5^Si canonical charmonium states. Although the mass of 
Y(4360) is somewhat smaller than the theoretical prediction, however, we notice that the 
mass predictions of various potential model for the high charmonium may differ by 10- 
lOOMeVSl], therefore Y(4360) as a 3'^Di cc state is not irrational. In this work we assume 
that the discrepancy in the spectrum is due to the theoretical uncertainties or other effects 
such as the coupled channel effects. It is interesting to refit the parameters a^, b, rric and a 
including both Y(4360) and Y(4660) or only Y(4660). 



III. ELECTRONMAGNETIC TRANSITIONS OF Y(4360) AND Y(4660) 



A. The e+e" leptonic decay of Y(4360) and Y(4660) 

The decay of quarkonium state into a lepton pair proceeds via a single virtual photon, 
as long as the mass of the initial quarkonim is sufficiently small that the contribution of a 
virtual Z can be ignored. The leptonic partial decay widths probe the compactness of the 
quarkonium system, and they provide useful information about the wavefunctions of the 
1 quarkonium states. The leptonic width of n'^Si charmonium is given byfl?!. [isl. 

r(.„%-.e-e-) = l^!^|p^(l-H^ + ^„) (5) 

where Cc = 2/3 is the charm quark electric charge, Mn is the mass of the n '^Si state, and the 
second term is the QCD correction. Vn(0) is the ri'^Si wavefunction at the origin, and the 
radial wavefunction is normalized according to Jq°° dr r^|'?/;„(r)p = 1. At the leading order, 
the width of D-wave cc states to e^e~ is proportional to 1-0^(0) p, 

r(n^D,-.e-e-) = ||||^:(0)r (6) 

which is generally smaller than the corresponding widths of the n'^Si states. Using the 
nonrelativistic quark model wavefunctions calculated in the previous section, we evaluate 
these leptonic decay widths of Y(4360) and Y(4660) at both the experimental values and the 
theoretical predictions of nonrelativistic potential model. The width predictions are given 



in Table [H where we choose ~ 0.23 12|. 



TABLE I: The e+e" partial widths of Y(4360) and Y(4660). 



Initial state 


Mass (GeV) 


r^+e- (keV) 


Y(4360)(33Di) 


4.361 


0.87 




4.455 


0.83 


Y(4660)(53Si) 


4.664 


1.34 




4.704 


1.32 



Using Eq.([T]), we can estimate that 

B(y(4360) Tx^Tx-^{2S)) ~ 1.20 x 10-^( or 1.36 x 10"^) 

B(y(4660) ^ 7r+7r-7/;(25)) ~ 2.24 x 10-^( or 5.67 x 10"^) (7) 
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The numbers in the bracket are the results corresponding to the solution II fit by the Belle 
collaboration. Generally we expect the branch fraction for cc(3^Di) — » 7r+7r~7/^(2S) should 
be of the order 10-^ e.g., i3(?/;(4160) 7r+7r-^(2^)) < 4 x IQ-^pj, where ?/^(4160) is a 
good candidate of 2^Di cc state. Therefore B{Y{4360) n'^7T~tlj{2S)) seems a little larger, 
which may be because of the QCD radiative corrections to r(Y(4360) — e^e~), and non- 
valence components may also contribute, which deserves investigating further. B{Y{A660) — >■ 
TT~^7r~ilj{2S)) ~ 10"'^, which indicates that Y(4660) may be a good candidate of 5 ^Si cc state. 



B. Radiative transitions of Y(4360) and Y(4660) 

Radiative decay of higher-mass charmonium states is an important way to produce lower 
charmonium states, and it plays significant role in charmonium physics. By means of the 
radiative transitions one can probe the internal charge structure of hadrons, hence it is 
useful for determining the quantum numbers and hadronic structures of heavy quark mesons. 
The radiative transition amplitude is determined by the matrix element of the EM current 
between the initial quarkonium state i and the final state /, i.e., (/li^mK)- Expanding in 
powers of photon momentum generates the electric and magnetic multipole moments, the 
leading order transition amplitudes are electric dipole (El) transion or magnetic dipole (Ml) 
transition. They are quite straightforward to be evaluated in the potential model. 



1. El transitions 



The partial width for El transitions between states n^^+^Lj and n'^^'+^Lj, cc state in 
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20 



21 



the nonrelativistic quark model is given by 
rEi(n^^+^L, ^n-^'^^L',, +7) = ^^(2J' + 1)5,. ^ K-"''^'LV|r|n ^^^^Lj) (8) 

where E^ is the photon energy, E, is the energy of final state n'^^'^^Lj,, and Mj is the 
mass of the initial state n^^+^Lj. We have included the relativistic phase factor and the 
statistical factor Sfi is 

{L' J' S I 
> (9) 
J L l/ 

The matrix element (n'^^ +^Lj, |r|n^^"*"^Lj) can be straightforwardly evaluated using the 
nonrelativistic Schrodinger wavefunctions of the model described in the previous section, 
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and the resulting El transition widths of Y(4360) and Y(4660) together with the photon 

energies are given in Table UTI- Table IVl where the El transition widths predictions for initial 

state assuming both the experiment observed masses and the nonrelativistic potential model 

predictions are given. The masses of the involved final state charmoniums are taken from the 
□ 

Particle Data Group [12| if the state is included in the meson summary table. If it is not, then 
the masses predicted in the nonrelativistic potential model described in the previous section 
are used. The exceptions are /ic(^Pi) and ?7c(3^So), we assume M(?7c(3 ^Sq)) = 4.011GeV(the 
mass of the known '?/;(4040) minus the theoretical 3S splitting) and M(/ic(^Pi)) = 3.525GeV, 
which is the spin-averaged mass of the ^Pj xj states. 

From Table UTI and Table UTTl we can see that Y(4360) should have very small El radiative 
widths to the triplet member of the IP multiplet and 1 state, if Y(4360) is a pure 3 '^Di 
state. The radiative widths to the unknow 3P and 2P triplet states xo(3'^Po), Xi(3'^Pi), 
Xo(2'^Po) and Xi(2'^Pi) are theoretically large, so the radiative decays of Y(4360) can be 
used to produce these states. Since the structures of both X(3940) and Y(3940) are still 
unclear, they possibly belong to the 2 '^Pj multiplet 0, 5]. Consequently the El transitions of 
Y(4360) into xo(2'^Po), Xi(2'^Pi) are especially of interests, which maybe helpful to clarifying 
the issue of X(3940) and Y(3940). 

Next we consider the El transition of Y(4660) as a 5 ^Si cc state. As is shown evidently in 
Table HVl and Table IVt the strong suppression of Y(4660) El decays to n'^Pj(n=l,2,3) states 
are predicted. The radiative width to 4'^Pj multiplet is large, which can provide access to 
the spin-triplet members of 4P multiplet. 

2. Ml transitions 

Ml transitions flip the quark spin, and Ml transitions are generally suppressed relative 
to the El transitions, and it has been observed in the charmonium system. Ml transition 
between different radial multiplets are only nonzero due to the small relativistic corrections 
to a vanishing lowest order Ml transition matrix element, therefore there may be serious 
inaccuracy in some channels. Analogous to the El transitions in the previous subsection. 



the Ml transitions width is given by 19|, l20|, l2l|, [22 1 



4r^F^E^ 9 T' -I- 1 V r V r 

(10)^ 



where the meaning of the notations is the same as that in the El transition case. The 
above formula has included the recoil factor jQ{E^r/2) with jo (2;) = sinx/x. Using the 
wavefunctions of nonrelativistic potential model in Sec. II, the Ml transitions width both 
with and without the recoil factor are calculated straightforwardly, theoretical predictions 
with the corresponding photon energies are shown in Table IVH -Table HXl Obviously the Ml 
transitions of Y(4360) and Y(4660) strongly depend on the recoil factors, and it may be to 
small to be observed. 



IV. STRONG DECAYS OF Y(4360) AND Y(4660) 

Strong decays of mesons are driven by nonperturbative gluodynamics, which are sensitive 
probe of hadron structure. However, it is very difficult to be calculated from the first 
principle. For charmonium above the DD threshold, the dominant decay modes usually 
are the open charm strong decays, in which the initial c and c separate into different final 
states. OZI forbidden decays are expected to be small, e.g. experimental indications are 
that i3(?^(3770) J/^ipTiTx) ~ 2.15 x 10"^ - 3.31 x 10-^ hence we shall focus on the open 
charm strong decays of Y(4360) and Y(4660) in this section. 

Although Open flavor decays are poorly understood from the QCD dynamics so far, a 
number of phenomenological models have been proposed to deal with this issue, the most 
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26|, the flux tube mo del |2 71. 



popular are the ^Pq model (^ark pair creation model) 23 

28] and the Cornell modeljl^, 0|- ■'■^ flux-tube model, a meson consists of a quark and 
antiquark connected by a tube of chromoelectric flux, which is treated as a vibrating string. 
For conventional mesons the string is in its vibrational ground state. The flux-tube breaking 



decay model 



28 



is similar to the ^Pq model, but extends it by including the dynamics of 
the flux tubes. This is done by considering the overlap of the flux tube of the initial meson 
with those of the two outgoing mesons. ^Pq model is a limiting case of the flux tube 
breaking model( ^Pq model emerges in the case of infinitely thick fiux tube ) , which greatly 



simplifies the calculations and gives similar results. The Cornell model [19|, |20|] assume that 
strong decays take place through pair production from the linear confining potential, which 
transform as the time component of a Lorentz vector rather than the Lorentz scalar 
in the ^Pq model. The Cornell model has the advantage of unifying the description of the 
spectrum and decays and completely specifies the strength of the decay. Recently it has 
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been used to discuss the possible charmonium assignments of X(3872) 29 1. 

The Orsay group have evaluated the open charm strong decays of three cc states ?/;(3770), 
'0(4040) and '?/;(4415) in the ^Pq model 30|], later these work was extended by taking into 



account flux tube breaking 



3l|. Recently T.Barnes et al used the ^Pq model to study the 



strong decays of both various candidates of X(3872) 32| and higher charmoinum up to the 
mass of the 4S multiplet 331] . In the following we shall consider the open flavor strong 
decays of Y(4360) and Y(4660) as 3^Di and 5^Si canonical charmonium in the simple 
harmonic oscillator wavefunction approximation in the framework of flux tube model, this 
approximation enables analytical studies of amplitudes, and it is known to be an excellent 
approximation for charmed mesons and light flavor mesons. Here we assume the harmonic 
oscillator parameter (3 of final states mesons are identical, which is different from (3a of the 
initial charmoinum. We will calculate the decay width following the procedure outlined in 



Ref. 



28 



34| . Previous attempts on exploring the charmonium strong decay in ^Pq model, flux 



tube model and Cornell model suggest that the typical error of the partial width predictions 
is 30%, and can reach factors of 2 or even 3. 

In the rest frame of A, the decay amplitude for an initial meson A into two final mesons 
B and C is, 



M{A^B + C) 



M 



(11) 



m + M' 

where both the flavor and spin overlap have been omitted in the above amplitude, and 
l{^Ai y) is the flux-tube overlap function, which measures the spatial dependence of the pair 
creation amplitude, y is the pair creation position, r^, and yc are respectively the quark- 
antiquark axes of A, B, and C mesons, they are related by = r^/2 + y, = r^/2 — y. 
The initial quark(antiquark) in A is of mass M with m the mass of the created quark pair. 
For charmonium decay concerned here, M = rric, m = mq(q=u,d,s). When the flux tube is 
in its ground states (conventional mesons), the flux-tube overlap function is 



28 



7irA,yj 



4° 

^00 




exp 



fb 2. 



(12) 



As usual, we take the string tension b = 0.18GeV^, and the constituent quark mass 
rriu = rrid = 0.33GeV, = 0.55GeV and = l.SGeV, and the estimated value / = 1.1 and 
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Ann = 1-0 are used in our calculation. The final D meson masses used to determined phase 

n 

space and final state momentum are taken from the Particle Data Group |12| and from recent 



Belle results 



351 ]. and if not available, the estimated mass motivated by the spectroscopy 



predictions are used(36|. These masses are M(D) = 1.8694GeV, M(D*) = 2.0078GeV, 
M(D*) = 2.308GeV(Belle), M(Di) = 2.444GeV, M(D;) = 2.422GeV, M(D2) = 2.459GeV, 
M(D') = 2.58GeV, M(D*') = 2.64GeV, M(Ds) = 1.9683GeV, M(D:) = 2.1121GeV, 
M(D*o) = 2.317GeV, M{D,i) = 2.459GeV, M(D^J = 2.535GeV, M(Ds2) = 2.572GeV, 
M(D^) = 2.67GeV, M(D*') = 2.73GeV. 

Heavy-light mesons are not charge conjugation eigenstates and so mixing can occur among 
states with the same J^. The = 1"*" axial vector cn and cs mesons Di and D'^^ are the 
coherent superpositions of quark model ^Pi and ^Pi states, 

|Di) = cos^l^Pi) +sin^|^Pi) 

\D[) = -sin^l^Pi) + cos^|=^Pi) (13) 

Little is known about the '^Pi- ^Pi mixing angle 6 at present, however, in the heavy quark 
limit, the mixing angle is predicted to be —54.7° or 35.3° if the expectation of heavy quark 



spin-orbit interaction is positive or negative 37 



. Since the former implies that the 2"*" 
state mass is larger than the 0+ state mass, and this agrees with experiment, we assume 
6 = —54.7° in the following. We note that generally finite quark mass will modify this 
mixing angle, and we can extract how large the mixing angle is by studying the dependence 
of the strong decay amplitudes on the mixing angle 6. 

When we calculate the decay widths from the amplitudes, there are ambiguities around 
the choice of phase space. The first choice is the fully relativistic phase space(RPS), which 
leads to a factor of in the final expression for the width in the center of mass frame, 
where Eb and Ec are respectively the energies of mesons B and C, and Ma is the mass of 
meson A. The second choice is fully non-relativistic phase space(NRPS), then the energy 
factor is replaced by ^^'^ , which is smaller than the relativistic phase space. A third 
possibility employed by Kokoski and Isgur, is the "mock meson" method, they suggest that 
the energy factor should be ^^sMc ^ ^^ere Mi(i=A,B,C) is the "mock meson" mass, which are 

Ma ' ' 



the calculated masses of the meson i in the spin-independent quark-antiquark potential 28 1. 
In practice, the numerical result is little different from the relativistic phase space except 
for the pseudoscalar goldstone bosons involved in the final states. Therefore we shall give 
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our partial width predictions for the relativistic phase space(RPS) and non-relativistic phase 
space(NRPS) in the following. 

Theoretical estimates for the harmonic oscillator parameters /? and /9a scatter in a relative 



26 



39 



and 



large region 0.3 — 0.7 GeV. Many recent quark model studies of mesons 
baryon 411] decays in ^Pq model use the value 0.4GeV. Moreover, the harmonic oscillator 
parameters of D, D* and D(^Pj) etc are predicted to be 0.45-0.66GeV, and mostly center 
around 0.5GeV|28|. Therefore we take /5a = 0.4GeV, (3 = O.SGeV in our calculation as an 
illustration, the outgoing mesons CM momentum p^, the partial widths and strong decay 
amplitudes for the kinematically allowed open charm decay modes of Y(4360) and Y(4660) 
are shown in Table |X] -Table IXVIIi We shall discuss some interesting and characteristic 
aspects about the strong decays of Y(4360) and Y(4660). 



A. Discussions about Y(4360) strong decay 

From Table IX]-Table IXIIl we can see that if Y(4360) is a pure 3^Di cc 
state with mass 4.361 GeV, there are ten open charm strong decay modes: 
DD,D*D,D*D*,D*D,D*D*,DiD,D'D,DsDs,D*Ds,D*D*, and the total width is predicted 
to be 67.69 MeV(RPS) or 53.24 MeV(NRPS), comparing with the Belle experimental mea- 
surement 74 ±15 ±10 MeV. Provided that Y(4360) mass is the prediction of non-relativistic 
potential model(4.455 GeV), then the additional decay modes DiD*, D'^D*, DsqD* and DsiDg 
become available. There is relative large difference between the relative phase space nor- 
malization and the non-relativistic phase space normalization in the DD mode, since the 
outgoing CM momentum ps is comparable to the D meson mass in this case. 

The leading decay mode of Y(4360) is predicted to be DD with a branching ratio~ 57%, 
the second-largest decay mode is D*D(^ 26%), and the D*D* mode also has sizable branching 
ratio. The relative partial wave amplitudes in the D*D* final state are predicted to have a 
very interesting pattern, = — -\/5, and Mf2 is predicted to be dominant, whereas it is 
zero for a S-wave charmonium decay. Measuring the relative branching ratio experimentally 
can determine whether Y(4360) is D-wave charmonium or S-wave charmoium. In addition, 
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we find the following relation, 

Msi(3 3Di ^ 'F,+ ^So) 

Msi(3=^Di ^ 3Pi+ '%) 

Therefore for the heavy mixing angle 6 

Msi(Y(4360) ^ Di + D) = Msi(Y(4360) ^ Di + D*) = (15) 

Thus, if Y(4360) is a pure 3=^01 cc state, the decays of Y(4360) to DiD or DiD*( if 
allowed by the phase space ) are in D-wave rather than in S-wave, where Di is the broader 
of the l"*" cq(q=u,d) axial vector mesons. To test the robustness of our conclusions, we 
study the stability of our results with respect to the variation of /3. The (3 dependence of 
the partial decay width and the total decay width are respectively shown in and FigH] and 
FigEl In FiglDwe showed the variation of DD, D*D, D*D*, D*D, D*D*, DiD and D'^D partial 
decay widths with the harmonic oscillation parameter j3, and we see that the partial decay 
widths into S+P final states are small. 



= i=Msi(3=^Di^ 'Pi+ 'Si) (14) 
= —54.7°, we have the following relations. 



_ DD 




d'd 




d'd' 









d'd' 




- D,D 




d'd 













0.40 0.42 0.44 0.46 0.48 0.50 0.52 0.54 0.56 0.58 0.60 

P(GeV) 




0.40 0.42 0.44 0.46 0.48 0.50 0.52 0.54 0.56 0.58 0.60 

p(GeV) 



FIG. 1: The variation of 

DD,D*D,D*D*,D*D,D*D*,DiD and D'-^D 
partial decay widths with /3 for Y(4360) as a 
3^Di charmonium state. 



FIG. 2: Y(4360) total width dependence on 

/3 as a 3 ^Di cc state in the flux tube model, 

the horizontal line denotes the current exper- 
□ 

imental upper bound [9[. 
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B. Discussions about Y(4660) strong decay 



Since the mass of Y(4660) is large, many open charm strong decay modes are allowable, 
which are listed obviously in Table IXTTTl -Table IXVIII For Y(4660) mass being the experimen- 
tal value 4.664 GeV, the total width is predicted to be 45.04 MeV(RPS) or 32.78 MeV(NRPS) 
by our parameters, which is in agreement with the Belle's measurement 48 ± 15 ± 3 MeV. 
Y(4660) dominantly decays into D*D, D*D* with branching ratios about 38% and 34% re- 
spectively, and DD is also an important mode. Some partial width of S+P final states can 
be comparable to the DD partial width for certain parameters. We would like to mention 
that Y(4660) as a 5 ^Si charmonium has four nodes in the radial wavefunction, consequently 
some modes of smaller branching ratios could be sensitive to the nodes positions. 

It is interesting to note that the flux tube model predicts the following relations between 
amplitudes, 

Msi(5% ^ 3Pi + ^Si) = -V2Msi(5 3Si ^ ^Pi + 

Mdi(5% ^ 3Pi + =^Si) = ^Mdi(5 3Si ^ iPi + ^Si) 

v2 

Md2(5=^Si ^ =^Pi + =^Si) = i=MD2(5^Si ^ iPi + =^Si) 

Msi(5 3Si ^ ig^) ^ -v^Mgi(5 3Si ^ iPi + ^So) 

Mdi(5%-^ 'Pi+ 'So) = -^Mdi(5 3Si^ ^Pi + ^So) (16) 

then the following interesting relation appears, 

Mdi(Y(4660) ^Di + D*) = — ^ Md2(Y(4660) ^ Di + D*) 

V3 

Mdi(Y(4660) ^ D; + D*) = — ^ Md2(Y(4660) ^ D', + D*) (17) 

v3 

The above two relations are independent of the ^Pi — ^Pi mixing angle 6. For the heavy 
quark mixing angle 6 = —54.7°, we have the following relations, 

Mdi(Y(4660) Di + D) = Mdi(Y(4660) Di + D*) = Md2(Y(4660) ^ Di + D*) = 
Msi(Y(4660) ^ D[ + D) = Msi(Y(4660) ^ D'^ + D*) = (18) 

The above relations imply that Y(4660) decays into both DiD and DiD* in S-wave, while 
into D^D and D'lD* in D -wave, if it is purely a 5 ^Si cs state. These predictions can also be 
used to test whether the ^Pi — ^Pi mixing is consistent with the heavy quark prediction. 
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As remarked in the previous discussion of Y(4360) decay, the D*D* mode is especially in- 
teresting. There are four partial wave amplitudes in this final state Mpo, Mpi, Mp2, Mf2, both 
Mpi and Mf2 are zero, and the ratio of two nonzero P-wave amplitude is Mp2/Mpo = — 2-\/5. 
However, this ratios is — 1/-\/5 in the case of D-wave charmonium decay, as is emphasized in 
the Y(4360) decay. Experimentally measuring these ratios are essential for understand the 
nature of Y(4660), and it is also a important test of the flux tube decay model. 

Moreover. Y(4660) can decay into D2D*, which is allowed by phase space. Five partial 
wave amplitudes are allowed for this process Msi, Mdi, Md2, M^s, Mqs, both Msi and Mqs 
amplitudes are predicted to be zero, whereas it is non-zero for a D-wave cc state decay. 
The ratios of the three D-wave amplitudes is Mm : Me)2 '■ M^s = 1 : '■ — 4-^/^. These 
predictions can be used to test whether Y(4660) is a S-wave charmonium, D-wave charmo- 
nium, or some other non-cc structure. In order to illustrate the parameter dependence of 
our predictions, we show the /3 dependence of the DD, D*D, D*D* and total S+P final states 
partial decay widths and the total width in Fig|3] and FigJH for Y(4660) as a 5 ^Si cc state. 
There are thirteen channels whose final states are S-wave and P-wave D mesons, and each 
partial decay width into S+P final state is at most close to the DD partial width for large 
part of the parameters regions. 
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FIG. 3: The variation of DD, D*D, D*D*and 
total S+P final states partial decay widths 
with P for Y(4660) as a charmonium 
state. 



FIG. 4: Y(4660) total width dependence on 
/3 as a 5 ^Si cc state in the flux tube model, 
the horizontal line denotes the current exper- 
imental upper bound[9[. 



In short summary, if Y(4360) and Y(4660) are 3 ^Di and 5 ^Si cc states respectively, the 
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DD, D*D and D*D* are expected to be the dominant decay modes, even if the variation of 
parameter is included. Some S+P final state may has comparable branching ratio for certain 
parameters. Careful measurements of these modes are crucial in testing these charmonium 
assignments. 



V. DISCUSSION AND CONCLUSION 



In this paper, we consider the canonical charmonium assignments of Y(4360) and Y(4660). 
Since these two structures are produced in initial state radiation from e"'"e~ annihilation and 
hence to have J^*" = 1 , so they are S-wave or D-wave states if they are canonical charmo- 
niums. From the mass spectrum prediction of the nonrelativistic potential model, we suggest 
Y(4360) is a 3^Di cc state, and y(4660) is a 5^Si cc state, if they are both conventional 
charmonium state, although Y(4360) mass is somewhat smaller than the nonrelativistic po- 
tential model prediction. We have investigated the e'^e~ leptonic decay. El transions. Ml 
transitions and open charm strong decay of both Y(4360) and Y(4660) in detail. 

Although the mass of Y(4360) is consistent with Y(4325) observed by the Babar col- 
laboration, it is much narrower. Thus more data are required to clarify whether they are 
the same structure. From the e^e~ partial width of Y(4360), we estimate B(Y {4360) — >■ 
7r+7r~?/'(25')) ~ 1.20 x 10~^( or 1.36 x 10~^), which is a little larger than the corresponding 
branching ratio for a conventional D-wave cc state decay. It possibly may be due to large 
QCD radiative corrections or some other non-cc componenets. It also indicates we should 



examine other possible interpretations of Y(^ 
effects especially deserve considering seriously 



360) further, the DiD and DsiDs threshold 



The Lattice QCD simulations predict that lightest charmonium h.jhnd{ccg) is about 



4.4GeV 



42 
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44| . It is obvious that Y(4360) mass is very near to 4.4GeV, so we can 



not completely exclude the possibility of Y(4360) as a 1 charmonium hybrid, a. 



though 



Y(4260) is already assumed to be a good candidate of 1 charmonium hybrid |45l. |46| . |47 |. 
Supposing that Y(4360) is a charmonium hybrid, its main decay modes should be 
D2D, DqD*, DiD, Dj^D according to the famous "S+P" selection rule in hybrid decay, and 
the DD,D*D,D*D* modes should be highly suppressed. Consequently measuring the 
DD, D*D, D*D*, D*D, D*D*, DiD, D^D modes are critical in distinguishing the canonical 
charmonium from the charmonium hybrid interpretation. 
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Xo(3^Po), Xi(3^Pi), Xo(2^Po) and Xi(2^Pi) are the main Y(4360) El transition modes 
as a 2^Di cc state, which possibly may be used to produce X(3940) and Y(3940), since they 
are expected to belong 2'^Pj multiplet. The strong suppression of Y(4660) El transitions to 
n^Pj(n=l,2,3) multiplet is predicted. The Ml transition of Y(4360) and Y(4660) should to 
too weak to be observed. 

We have discussed the open charm strong decays of Y(4360) and Y(4660) in the flux 
tube model in detail. Both Y(4360) and Y(4660) are predicted to dominantly decay into 
DD, D*D, D*D*, the partial width of some S+P final states can be comparable with that 
of DD, D*D or D*D* for certain parameters. Measuring the ratios of the amplitudes in the 
D*D* final state will show whether Y(4360) and Y(4660) are consistent with the charmonium 
assignments made in this work. If Y(4360) is a pure 3 ^Di charmonium, Mf2 amplitude is 
predicted to be largest and = For Y(4660) as a pure 5'^Si cc state, we predict 

that the amplitude Mf2 is zero and = —2\/5. Similarly the D2D* amplitude ratios 
in Y(4660) decay can test whether Y(4660) is a S-wave or D-wave cc state, although the 
branching ratio of D2D* is predicted to be small. Provided that Y(4660) is a 5 ^Si cc state, 
we have Mdi : Md2 : Mds = 1 : y'^ : — 4y^, and the amplitudes Msi and Mgs is zero, 
which is non-zero for a D-wave state decay. The above results are generally correct for 
S-wave or D-wave initial state decay. The careful measurement of these relative branching 
ratios would play a critical role in understanding Y(4360) and Y(4660). 

The Belle and Babar Collaboration have measured the exclusive e"'"e~ — DD, e^e~ 
DD* and e"'"e~ D*D* cross section using initial state radiation 
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50, 



511 ] ■ and the shapes 



of the cross sections are similar. There is a peak in the Y(4660) region, however, no structure 
is clearly observed near position of Y(4360) so far. Therefore Y(4660) as a 5^Si cc state 
is consistent with current experimental data, however, Y(4360) may be a state beyond the 
quark model. Since the DDi threshold is 4291.4 MeV (Md + Md^ 4291.4MeV), which 
is close to Y(4260), a possible way of reconciling Y(4260) and Y(4360) is that Y(4260) is 
mainly the S-wave DDi threshold effect and Y(4360) is a charmonium hybrid. The relevant 
work is in progress 52l| . 

The confirmations and more experimental studies on Y(4360) and Y(4660) at BES and 
CLEO are expected. Careful study of Y(4360) and Y(4660) will greatly shed light on the 
charmonium spectroscopy. 



16 



ACKNOWLEDGEMENTS 

We acknowledge Prof. Dao-Neng Gao for very helpful and stimulating discussions, and 
we are especially grateful to Prof. Chang-Zheng Yuan for his professional comments. This 
work is partially supported by National Natural Science Foundation of China under Grant 
Numbers 90403021, 10005008, and KJCX2-SW-N10 of the Chinese Academy. 



[1] J. L. Rosner et al.[CLEO Collaboration], Phys. Rev. Lett. 95, 102003 (2005) , hep-ex/0505073 



[2] S. K. Choi et a/. [BELLE collaboration], Phys. Rev. Lett. 89, 102001 (2002) [Erratum-ibid. 



89, 129901 (2002)], [hep-ex/0206002| P. Rubin et a/.[CLEO Collaboration], Phys. Rev. D 72, 



092004 (2005), hep-ex/0508037 



[3] S. K. Choi et aZ.[Belle Collaboration], Phys. Rev. Lett. 91, 262001 (2003), hep-ex/0309032 



D. Acosta et a/.[CDF II Collaboration], Phys. Rev. Lett. 93, 072001 (2004), hep-ex/0312021 



V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 93, 162002 (2004), hep-ex/0405004 



B. Aubert et aL[BABAR Collaboration], Phys. Rev. D 71, 071103 (2005), hep-ex/0406022 



K. Abe et a/. [Belle Collaboration], hep-ex/0507019 



K. Abe et a/. [Belle Collaboration], Phys. Rev. Lett. 94, 182002 (2005), hep-ex/0408126 



K. Abe et a/.[Belle Collaboration], |hep-ex/0507033 S. Uehara et aZ.[Belle Collaboration], 



Phys. Rev. Lett. 96, 082003 (2006), hep-ex/0512035 



B. Aubert et al.[B ABAR Collaboration], Phys. Rev. Lett. 95, 142001 (2005), hep-ex/0506081 



E. S. Swanson, Phys. Rept. 429, 243 (2006), |hep-ph/060lTl()| E. Eichten, S. Godfrey 

H. Mahlke and J. L. Rosner, |hep-ph/0701208| 

X. L. Wang et aZ.[Belle Collaboration], |ar2£iYi070L36991hep-ex] . 



B. Aubert et aZ.[BABAR Collaboration], Phys. Rev. Lett. 98, 212001 (2007), [hep^ex/0610057 

K. Abe et a/. [The Belle Cohaboration], |M2LiYi0708T790J [hep-ex]. 

Particle Data Group, W.-M. Yao et al, Journal of Physics G 33, 1 (2006). 

N. Brambilla, A. Pineda, J. Soto and A. Vairo, Rev. Mod. Phys. 77, 1423 (2005), 

arXiv:hep-ph/041 0047[ 



N. Brambilla and A. Vairo, |arXiv:hep-ph/9904330 



T. Barnes, S. Godfrey and E. S. Swanson, Phys. Rev. D 72, 054026 (2005) , hep-ph/0505002 



17 



[16 

[17: 

[18 
[19 

[20 

[21 

[22; 

[23 
[24 

[25; 
[26 

[27; 

[28 
[29 
[30 

[31 

[32; 

[33 
[34 
[35; 
[36 

[37; 

[38 
[39 



W. Lucha and F. F. Schoberl, Int. J. Mod. Phys. C 10, 607 (1999), |hep-ph/98H453 



R. Van Royen and V. F. Weisskopf, Nuovo Cim. A 50 (1967) 617[Erratum-ibid. A 51 (1967) 
583]. 

W. Kwong, J. L. Rosner and C. Quigg, Ann. Rev. Nucl. Part. Sci. 37, 325 (1987). 

E. Eichten, K. Gottfried, T. Kinoshita, J. B. Kogut, K. D. Lane and T. M. Yan, Phys. Rev. 

Lett. 34, 369 (1975) [Erratum-ibid. 36, 1276 (1976)]. 

E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane and T. M. Yan, Phys.Rev.Lett.36:500,1976; 
Phys. Rev. D 17, 3090 (1978) [Erratum-ibid. D 21, 313 (1980)]; Phys.Rev.D21:203,1980. 
W. Kwong and J. L. Rosner, Phys. Rev. D 38, 279 (1988). 



N. Brambilla et al. [Quarkonium Working Group], arXiv:hep-ph/0412158 
L. Micu, Nucl. Phys. B 10, 521 (1969). 

A. Le Yaouanc, L. OUver, O. Pene and J. C. Raynal,Phys. Rev. D 8 (1973) 2223; Phys. Rev. 
D 9, 1415 (1974); Phys.Rev.Dll:680,(1975); Phys. Rev. D 11, 1272 (1975). 



P. Geiger and E. S. Swanson, Phys. Rev. D 50, 6855 (1994), hep-ph/9405238 



E. S. Ackleh, T. Barnes and E. S. Swanson, Phys. Rev. D 54, 6811 (1996), |hei>ph/9604355 
N. Isgur and J. E. Paton, Phys. Rev. D 31, 2910 (1985). 
R. Kokoski and N. Isgur, Phys. Rev. D 35, 907 (1987). 



E. J. Eichten, K. Lane and C. Quigg, Phys. Rev. D 69, 094019 (2004), [hep-ph/0401210 

A. Le Yaouanc, L. Ohver, O. Pene and J. C. Raynal, Phys. Lett. B 71, 397 

(1977);Phys.Lett.B72:57,1977. 



P. R. Page, Nucl. Phys. B 446, 189 (1995), hep-ph/9502204 



T. Barnes and S. Godfrey, Phys. Rev. D 69, 054008 (2004), [hep^ph/031 1162 



T. Barnes, S. Godfrey and E. S. Swanson, Phys. Rev. D 72, 054026 (2005), hep-ph/0505002 



G. J. Ding and M. L. Yan, Phys. Lett. B 650, 390 (2007), hep-ph/0611319 ; hep-ph/0701047 



K. Abe et al. [Belle Collaboration], Phys. Rev. D 69, 112002 (2004), [hei>ex/0307021 

S. Godfrey and N. Isgur, Phys. Rev. D 32 (1985) 189. 

S. Godfrey and R. Kokoski, Phys. Rev. D 43, 1679 (1991). 



S. Godfrey, Phys. Rev. D 72, 054029 (2005), hep-ph/0508078 



T. Barnes, F. E. Close, P. R. Page and E. S. Swanson, Phys. Rev. D 55, 4157 (1997), 



|hep-ph/9609339 



[40] T. Barnes, N. Black and P. R. Page, Phys. Rev. D 68, 054014 (2003), nucl-th/0208072 



18 



S. Capstick and W. Roberts, Phys. Rev. D 49, 4570 (1994), |nucl-th/9310030 



C. W. Bernard et al. [MILC Collaboration], Phys. Rev. D 56, 7039 (1997), hep-lat/9707008 



Z. H. Mei and X. Q. Luo, Int. J. Mod. Phys. A 18, 5713 (2003), jhep-lat/0206012 



G. S. Bah, Eur. Phys. J. A 19, 1 (2004), hep-lat/0308015 



F. E. Close and P. R. Page, Phys. Lett. B 628, 215 (2005), hep-ph/0507199 



S. L. Zhu, Phys. Lett. B 625, 212 (2005), hep-ph/0507025 



E. Kou and O. Pene, Phys. Lett. B 631, 164 (2005), [hep^h/0507119 
G. J. Ding. larXiv: 071 1.1485 1 [hep-oh]. 



K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 98, 092001 (2007), |hep-ex/0608018 
G. Pakhlova, larXiv:0708.0082] [hep-ex]. 

T. B. Collaboration and B. Aubert. [arXiv:0710.137T] [hep-ex]. 
Gui-Jun Ding and Mu-Lin Yan, work in progress. 



19 



TABLE II: El radiative transitions of Y(4360)(3^Di), and the Y(4360) mass is taken from exper- 
iment. 

Final meson E^ (MeV) T (keV) 

X2(33P2) 44 0.12 

Xi(3^Pi) 89 15.9 

Xo(33Po) 156 112 

X2(2 3P2) 371 0.39 

Xi(2^Pi) 414 8.00 

Xo(2=^Po) 479 16.2 

X2(2^F2) 10 0.062 

X2(1^P2) 731 0.17 

Xi(l^Pi) 767 2.94 

843 5.11 



Xo(l=^Po) 

X2(1^F2) 319 



0.053 
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TABLE III: El radiative transitions of Y(4360)(33Di), and the Y(4360) mass is the prediction of 



the nonrelativistic potential model, which is 4.455GeV. 

Final meson E^ (MeV) F (keV) 

X2(3^P2) 135 3.65 

Xi(33Pi) 180 128 

Xo(3 3Po) 246 425 

X2(2 3P2) 456 0.71 

Xi(23Pi) 498 13.7 

Xo(2 3Po) 562 25.6 

X2(2^F2) 103 58.3 

X2(1^P2) 808 0.23 

Xi(l^Pi) 844 3.85 

Xo(l^Po) 918 6.49 

X2(1^F2) 405 0.11 



21 



TABLE IV: El radiative transitions of Y(4660)(53Si), and the Y(4360) mass is taken from exper- 
iment. 



Final meson E^ (MeV) T (keV) 

X2(43P2) 42 7.92 

Xi(43Pi) 87 42.8 

Xo(43Po) 152 75.5 

X2(3 3P2) 334 0.34 

Xi(3 3Pi) 377 0.29 

Xo(3=^Po) 439 0.15 

X2(2 3P2) 640 0.68 

Xi(2 3Pi) 681 0.48 

Xo(2=^Po) 741 0.20 

X2(1^P2) 976 0.37 

Xi(l^Pi) 1010 0.24 

Xo(l^Po) 1082 0.097 
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TABLE V: El radiative 


transitions of Y(4660)(5^Si), and the Y(4660) mass 


is the prediction of 


the nonrelativistic potential model, which 


is 4.704GeV. 






H m 51 1 TTi ciG/~m 




E^ (MeV) 




r (keV) 






ol 




57.7 


Xi(4=^Pi) 








129 


Xo(4 3Po) 




191 




147 


x^^f'i 3p \ 




oil 




0.46 


Xi(3=^Pi) 








0.38 


Xo(33Po) 




475 




0.19 






d75 




0.78 






714 




0.55 


Xo(2=^Po) 




774 




0.23 


X2(1=^P2) 




1008 




0.40 


Xi(l=^Pi) 




1042 




0.26 


Xo(l=^Po) 




1112 




0.10 


TABLE VL Ml radiative transitions of Y(4360)(3 ^Di) 


, and Y(4360) mass 


is the experimental 


value 4.361GeV. 










Final meson 


E^ (MeV) 




r (keV) 


Tree (keV) 


??c2(2^D2) 


199 




0.00031 


0.082 


??c2(l'D2) 


525 




0.00067 


0.19 


TABLE VII: Ml radiative transitions of Y(4360)(3 ^Di), 


and Y(4360) mass is 


the theoretical pre- 


diction of the nonrelativistic potential mo 


del, which is 4.455GeV. 




Final meson 


E^ (MeV) 




r (keV) 


Trec (keV) 


?/c2(2^D2) 


287 




0.00093 


0.24 


r?c2(l'D2) 


607 




0.0010 


0.29 
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TABLE VIII: Ml radiative transitions of Y(4660)(5 ^Si), and Y(4660) mass is the experimental 
value 4.664GeV. 



Final meson (MeV) T (keV) Tree (keV) 

7/c(4iSo) 272 0.15 0.95 

7?c(3iSo) 607 0.47 3.45 

7?c(2iSo) 913 0.82 4.26 

7?c(l^So) 1381 2.65 9.36 



TABLE IX: Ml radiative transitions of Y(4660)(5 ^Si), and Y(4660) mass is the theoretical pre- 
diction of the nonrelativistic potential model, which is 4.704GeV 

Final meson E^ (MeV) T (keV) Tree (keV) 

r/c(5^So) 19 0.013 0.013 

7/c(4^So) 309 0.23 1.39 

7?c(3^So) 642 0.55 4.05 

7?c(2^So) 945 0.90 4.69 

7?c(l^So) 1409 2.80 9.89 
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TABLE X: Open-charm strong decay of Y(4360)(3'^Di), a factor of +i has been suppressed in all 
old partial waves. Y(4360) mass is the experimental value 4.361GeV. 



Mode PB (GeV) T (MeV) Amps(GeV-V2) 







RPS 


NRPS 




DD 


1.12 


39.03 


28.69 


Mpo = 0.3877 


D*D 


1.00 


18.02 


14.24 


Mpi = -0.1977 


D*D* 


0.85 


6.83 


5.79 


Mpo = 0.0705 
Mpi = 
Mp2 = -0.0315 
Mf2 = 0.1695 


DP 


0.26 


0.26 


0.25 


Md2 = -0.0463 


D*D* 


0.31 


0.67 


0.66 


Msi =0 

Mdi = -0.0684 


DiD 


0.32 


0.62 


0.60 


Msi =0 

Mdi = -0.0653 


D'iD 


0.38 


1.56 


1.50 


Msi = 0.0776 
Mdi = -0.0534 


DsDs 


0.94 


1.23 


1.01 


Mpo = 0.1066 


D*Ds 


0.77 


0.00 


0.00 


Mpi = 4.3 X 10- 


DsDs 


0.54 


0.11 


0.10 


Mpo = -0.0135 
Mpi = 
Mp2 = 0.0060 
Mf2 = -0.0384 



Total 68.33 52.84 
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TABLE XI: Open-charm strong decay of Y(4360)(3^Di), Y(4360) mass is 4.455GeV the prediction 
of the nonrelativistic potential model. 



Mode 


PB (GeV) 


RPS 


r (MeV) 

NRPS 


Amps(GeV-i/2) 


DD 


1.21 


51.33 


36.15 


Mpo = 0.4234 


D*D 


1.10 


35.42 


26.81 


Mpi = -0.2615 


D*D* 


0.96 


18.22 


14.80 


Mpo = 0.1411 










Mpi = 










Mp2 = -0.0631 










Mf2 = 0.2367 




0.52 


0.02 


0.02 


Md2 = -0.0089 




0.55 


0.00 


0.00 


Msi =0 










Mdi = 0.0034 


DiD* 


0.08 


0.01 


0.01 


Msi =0 










Mdi = 0.0117 










Md2 = 0.0003 


DiD 


0.55 


0.00 


0.00 


Msi =0 










Mdi = 0.0004 




0.24 


0.50 


0.49 


Msi = -0.0323 










Mdi = -0.0304 










Md2 = -0.0506 




0.59 


0.34 


0.31 


Msi = -0.0046 










Mdi = 0.0345 


DsDs 


1.04 


4.70 


3.67 


Mpo = 0.1953 


DsDs 


0.89 


0.55 


0.46 


Mpi = -0.0509 


DsDs 


0.71 


0.02 


0.02 


Mpo = -0.0117 










Mpi = 










Mp2 = 0.0052 










Mf2 = 0.0091 
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TABLE XII: Open-charm strong decay of Y(4360)(3^Di), Y(4360) mass is 4.455GeV the prediction 
of the nonrelativistic potential model(continued). 



Mode 


PB (GeV) 


RPS 


r (MeV) 

NRPS 


Amps(GeV-i/2) 




0.24 


0.12 


0.12 


Msi =0 










Mdi = -0.0462 


DsiDs 


0.25 


0.11 


0.10 


Msi =0 



Mdi = -0.0428 

Total 111.35 82.96 
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TABLE XIII: Open-charm strong decay of Y(4660)(5^Si). Mlj is the partial wave amplitude, 
where L=S, P, D,... is the relative angular momentum and J is their total spin. Note that a 
factor of +i has been suppressed in all old partial waves. Y(4660) mass is the experimental value 

4.664GeV. 

Mode pB (GeV) T (MeV) Amps(GeV"i/2) 







RPS 


NRPS 




DD 


1.39 


5.29 


3.40 


Mpo = 0.1238 


D*D 


1.30 


17.47 


12.06 


Mpi = 0.1651 


D*D* 


1.19 


15.32 


11.36 


Mpo = 0.0499 
Mpi = 
Mp2 = -0.2230 
Mf2 = 


D*D* 


0.67 


0.23 


0.21 


Msi =0 
Mdi = 0.0042 
Md2 = 0.0055 
Md3 = -0.0258 
Mg3 = 


D*D 


0.86 


0.75 


0.64 


Md2 = -0.0422 




0.88 


1.00 


0.86 


Msi = -0.0480 
Mdi = 


DiD* 


0.69 


0.22 


0.20 


Msi = -0.0256 
Mdi = 
Md2 = 


DiD 


0.88 


0.94 


0.80 


Msi = -0.0467 
Mdi = 


b[d* 


0.73 


0.01 


0.01 


Msi =0 

Mdi = -0.0032 

Md2 = 0.0056 




0.91 


1.78 


1.50 


Msi =0 
Mdi = 0.0634 


d*'d* 


0.19 


0.54 


0.54 


Mpo = 0.0166 
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TABLE XIV: Open-charm strong decay of Y(4660)(53Si). Y(4660) mass is the experimental value 
4.664GeV(continued). 



MnHe 


r)D (TleVl 




r CMeVl 


ATTms^^;pv~V2^ 

iT-lll^Dl VjC V J 


n*'n* 








Mr,M — 










Mt^o — —0 074 

iVip2 — U.U('4:0 










i\/T__ — n 


T)*'T) 


\j .Oct 




9fi 


Mm — —0 031 3 


DD* 


0.42 


n Dfi 


05 


Mdi — 01 fi5 


d'd 




n 04 


04 


Mnn — 01 07 

1V±J-'Q U.UXU 1 


D D 


1.25 


0.50 


0.35 


Mr>n — 0567 






"^Q 


30 


Mr>i — 0373 

IVJ.J-']^ — U.UO 1 o 




QQ 


no 

u.uu 


DO 
u.uu 


iVipQ — U.UUX'4: 










Moi — 










Mdo — —0 0061 

IV J. J-* 2 — u . u u u ± 










Mtpo — 




VJ . (JO 


n on 


DO 


Mt^o — 0031 


^sO^s 




Ofi 


05 

U.UU 


Moi — 01 80 

1V±J;-j2 — u . u ± ou 










Mm — 


J-^slJ-'s 


n 4fi 


09 


09 


Moi — 01 91 










MrM — 










Mm — 






Ofi 


05 

U.UU 


Moi — 01 80 










MrM — 


^sl^s 




01 

U.U 1 


01 

U.U-L 












Mdi = 0.0052 










Md2 = -0.0090 


D^iDs 


0.60 


0.02 


0.02 


Msi =0 










Mdi = -0.0106 


dIDs 


0.24 


0.05 


0.05 


Mpo = 0.0280 


Total 




45.04 


32.78 
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TABLE XV: Open-charm strong decay of Y(4660)(5^Si). Mlj is the partial wave amphtude, where 
L=S, P, D,... is the relative angular momentum and J is their total spin. Note that a factor of +i 
has been suppressed in all old partial waves. Y(4660) mass is the prediction of the nonrelativistic 
potential model, which is 4.704GeV. 



Mode 


Or (GeV) 


RPS 


r (MeV) 

NRPS 


Amps(GeV-V^) 


DD 


1.43 


5.13 


3.24 


Mpn = 1200 


D*D 


1.33 


19.82 


13.45 


Mpi = 1728 


D*D* 


1.23 


20.93 


15.25 


Mpn = 0571 










Mpi = 










Mpo — —0 2554 










Mwo — 




0.73 


0.01 


0.01 


Mqi = 










Mm — 0009 










Mn9 = 0012 










Md3 = -0.0057 












DSD 


0.91 


2.00 


1.68 


Mn9 = — 0666 


DSD* 


0.93 


3.37 


2.83 


Mci — —0 0852 










Mm = 


DiD* 


0.76 


0.06 


0.05 


Mci = -0 0122 










Mdi = 










Md2 = 


DiD 


0.93 


3.19 


2.67 


Msi = -0.0834 










Mdi = 




0.79 


0.12 


0.11 


Msi =0 










Mdi = 0.0089 










Md2 = -0.0155 


d'^d 


0.96 


3.54 


2.93 


Msi =0 



Mdi = 0.0865 



30 



TABLE XVI: Open-charm strong decay of Y(4660)(52Si), Y(4660) mass is 4.704GeV the prediction 
of the nonrelativistic potential model(continued). 



Mode 






r fMeV) 


AmDsfGeV"^/^) 


r)*'n* 


u.ou 


01 
u.ux 


01 
u.ux 


iVlpQ — U.uUXO 










Mr>i — n 










Mr>n — — n nn5s 

ivip2 — u.uutjo 










ivxj-< 2 — u 




u.uu 


DO 
u.uu 


DO 
u.uu 




DD* 








Mdi — 0'^79 


d'd 




n fiQ 


fi9 


Mor, — 04SS 


D D 


1.29 


0.70 


0.49 


Mr>n — 0659 




1 1 7 


n 7fi 

U. 1 u 


57 


Mm — 0509 

iVXp2 U.UUUt? 




-L . WO 


u.uu 


05 

U.UU 


iVXpQ U. UU'iO 










Mdi — 










Mnn — —0 091 7 

lVXj--*2 U.UZ;X { 










Mt?o — 


V)* V)* 




u.ux 


01 

u.ux 


Mcl — 










MrM — —0 009fi 

1VXJJ2 U.UUiiiU 










Mno — —0 00S4 










iVT — n ni fin 

ivip)3 — u.uiuu 














n fii 


n 09 


09 


Mt-,o — 01 1 


d:oDs 


0.79 


0.05 


0.04 


Msi = 0.0159 










Mdi =0 


DsiD* 


0.55 


0.00 


0.00 


Msi = -4.3 X 10-5 










Mdi =0 










Md2 = 


DsiDs 


0.79 


0.05 


0.04 


Msi = 0.0160 



Mdi =0 
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TABLE XVII: 


Open-charm strong 


decay of Y(4660)(53Si), Y(4660) 


mass IS 4.7U4GeV the predic- 


tion of the nonrelativistic potential 


model(continued) . 






l\/rnrlo 


Pb ['^ev ) 


i (^iViev ) 




Amps(GeV-i/2) 






tiro 








n "if; 
U.oD 


n no 
U.Uz 


n no 
U.Uz 


Msi =0 










Mdi = 0.0073 










Md2 = -0.0126 


DliDs 


0.67 


0.04 


0.03 


Msi =0 










A /r n ni /I n 
iVL£)i = — U.Ui4y 


Ds'Ds 


0.11 


0.03 


0.03 


Mpi = 0.0314 


DlD, 


0.39 


0.00 


0.00 


Mpo = 0.0051 


Total 




60.96 


44.48 
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